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ABSTRACT 

In  a  partially  ionized  hypersonic  flow  field,  electron  collision  plays  a  significant  role  in  energy  re-distribution 
and  in  the  production  of  excited  states  of  atoms  and  molecules  that  can  subsequently  radiate.  Thus  electron 
collision  data  is  part  of  the  database  used  in  modelling  high  speed  entry  flows.  First  principles  computational 
methods  for  the  calculation  of  electron-impact  data  are  reviewed.  Their  merits  are  discussed,  based  on  the 
requirements  of  reliability  and  the  ability  to  handle  high-lying  excited  states.  Three  processes  are  considered: 
the  production  of  electrons  by  electron-impact  ionization,  the  removal  of  electrons  by  radiative  recombination 
and  dielectronic  recombination,  and  the  production  of  electronic  excited  states  by  electron  impact.  Electron  - 
nitrogen  atom  collisions  are  used  as  an  illustrative  example.  The  study  clarifies  the  role  of  high-lying 
electronic  excited  states  in  the  ionization  process.  It  shows  that  the  recombination  processes  can  produce 
sizable  radiative  heat  load  in  the  VUV  region.  It  also  shows  examples  using  a  scaling  method  to  improve  the 
frequently  used  Born  approximation  for  electron  -impact  excitations. 


1.0  INTRODUCTION 

During  the  hypersonic  entry  of  a  space  vehicle  into  a  planetary  or  lunar  atmosphere,  the  flow  field  becomes 
partially  ionized.  The  percentage  of  ionization  depends  on  the  entry  speed  and  the  vehicle  size.  In  this  regime, 
electron  collision  provides  an  efficient  means  of  producing  electronic  excited  states  of  the  atoms  and 
molecules  in  the  flow  field.  Thus  it  plays  an  important  role  in  determining  the  internal  energy  and  state 
distribution  of  the  gaseous  particles.  The  excited  states  in  turn  are  the  source  of  radiation  observed  during  a 
hypersonic  entry. 

Electron-atom/molecule  collisions  differ  from  heavy  particle  (atom-atom,  atom-molecule  and  molecule- 
molecule)  collisions  in  two  aspects.  First,  the  mass  of  an  electron  is  more  than  four  orders  of  magnitude 
smaller  than  the  reduced  mass  of  N2.  Thus  its  average  speed,  and  hence  its  average  collision  frequency,  is 
more  than  100  times  larger.  Even  in  the  slightly  ionized  regime  with  only  1%  electrons,  the  frequency  of 
electron-atom/molecule  collisions  is  equal  to  or  larger  than  that  of  heavy  particle  collisions,  an  important 
consideration  in  the  low  density  part  of  the  atmosphere  where  the  reaction  probability  is  frequently  controlled 
by  the  collision  frequency.  Second,  the  interaction  potential  between  a  charged  particle  (electron)  and  a 
neutral  particle  is  longer  range  than  neutral-neutral  interactions.  Hence  electron-atom/molecule  collision  cross 
sections  tend  to  be  larger.  Also,  low-energy  electron  collisions  can  be  effective  in  spin  changing  excitations. 
Another  characteristic  of  electron  collisions  is  that  it  generally  produces  a  variety  of  excited  states  whereas 
heavy  particle  collisions  tend  to  produce  specific  excited  states. 

Modelling  electron  collisions  in  nonequilibrium  gas  dynamics  requires  data  to  simulate  its  production  and 
removal  in  the  flow  field.  Also,  in  the  non-equilibrium  regime  where  the  populations  of  radiative  species  are 


RTO-EN-AVT-162 


6-1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

SEP  2009 


2.  REPORT  TYPE 

N/A 


3.  DATES  COVERED 


4.  TITLE  AND  SUBTITLE 

Electron-Impact  Excitation  and  Ionization  in  Air 


6.  AUTHOR(S) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 
5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Huo  Consulting  LLC,  Los  Altos,  CA  94024-3168  USA  REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADA562449.  RTO-EN-AVT-162,  Non-Equilibrium  Gas  Dynamics  -  From  Physical  Models  to 
Hypersonic  Flights  (Dynamique  des  gaz  non-  equilibres  -  Des  modeles  physiques  jusqu’au  vol 
hypersonique).,  The  original  document  contains  color  images. 

14.  ABSTRACT 

In  a  partially  ionized  hypersonic  flow  field,  electron  collision  plays  a  significant  role  in  energy 
re-distribution  and  in  the  production  of  excited  states  of  atoms  and  molecules  that  can  subsequently 
radiate.  Thus  electron  collision  data  is  part  of  the  database  used  in  modelling  high  speed  entry  flows.  First 
principles  computational  methods  for  the  calculation  of  electron-impact  data  are  reviewed.  Their  merits 
are  discussed,  based  on  the  requirements  of  reliability  and  the  ability  to  handle  high-lying  excited  states. 
Three  processes  are  considered:  the  production  of  electrons  by  electron-impact  ionization,  the  removal  of 
electrons  by  radiative  recombination  and  dielectronic  recombination,  and  the  production  of  electronic 
excited  states  by  electron  impact.  Electron-nitrogen  atom  collisions  are  used  as  an  illustrative  example.  The 
study  clarifies  the  role  of  high-lying  electronic  excited  states  in  the  ionization  process.  It  shows  that  the 
recombination  processes  can  produce  sizable  radiative  heat  load  in  the  VUV  region.  It  also  shows  examples 
using  a  scaling  method  to  improve  the  frequently  used  Born  approximation  for  electron-impact  excitations. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

SAR 

22 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Electron-Impact  Excitation  and  Ionization  in  Air 


e-t&r-iA-^rt'v!  Artt/rt  kt 


determined  by  the  quasi-steady-state  approximation,  the  electron  number  density,  temperature,  and  electron- 
impact  excitation  rate  coefficients  are  part  of  the  input  data  for  QSS  calculation.  Similar  data  needs  also  exist 
in  other  disciplines  dealing  with  plasma,  including  astrophysics,  fusion  science,  plasma  processing  in 
microelectronics,  etc.  While  the  composition  and  temperature  of  the  plasmas  may  be  different,  these 
communities  frequently  employ  a  common  set  of  approximate  formulas  for  the  cross  section/rate  coefficient 
expressions.  Using  a  combination  of  experimental  data  and  a  variety  of  approximate  formulas,  a  number  of 
databases  have  been  developed  in  the  entry  physics  community.  The  data  sets  by  Park  [1,2],  Losev  [3]  and 
Bird’s  TCE  [4,5]  models  are  well  established  examples  of  this  approach.  The  NEQAIR  package  [6]  that 
simulates  nonequilibrium  radiation  in  an  entry  flow  employs  Gryzinski’s  [7]  classical  formula  is  used  for 
electron-impact  excitation  of  atoms  whereas  for  molecules  experimental  data  are  used,  sometimes  by 
extrapolation  or  by  analogy.  More  recent  models  have  incorporated  improved  databases,  based  on  new 
experimental  data  and/or  theoretical  calculations.  The  collisional  radiative  model  by  Bourdon  et  al.  [8] 
includes  many  updates.  Similarly  SPRADIAN07  [9]  incorporated  new,  improved  data  into  the  NEQAIR 
model. 

As  part  of  a  systematic  development  of  a  physics-based  model  of  nonequilibrium  chemistry  and  radiation  in 
hypersonic  flow,  it  is  timely  to  investigate  and  update  the  electron  collision  cross-sections/rate  coefficients 
based  on  first  principles  calculations.  The  present  focus  is  on  electron-atom  collisions,  because  atomic 
radiation  dominates  at  high  entry  speed.  Our  study  covers  electron-impact  ionization.  Together  with 
associative  ionization,  they  are  the  production  mechanisms  of  electrons  in  the  flow  field.  Two  electron 
removal  mechanisms  are  investigated,  radiative  recombination  and  dielectronic  recombination.  Selected 
electron-impact  excitations  of  bound-bound  transitions  are  also  studied. 


2.0  QUANTUM  MECHANICAL  FORMULATION  OF  ELECTRON  COLLISIONS 

Due  to  the  light  electron  mass,  in  general  e-atom/molecule  collisions  should  be  treated  using  quantum 
mechanics  instead  of  classical  mechanics.  The  Schrodinger  equation  for  an  e-atom  system  is  given  by 

(//-£TF(r,L  rv+1)  =  0  (1) 

Here  H  is  the  Hamiltonian  of  the  e  +  atom  system,  E  the  total  energy,  T  the  corresponding  wave  function,  and 
Ti  the  spatial  and  spin  coordinates  of  the  itb  electron.  The  indices  i  and  j  label  the  bound  electrons,  and  N+l 
label  the  free  electron.  Due  the  large  difference  between  electron  and  nuclear  mass,  it  is  assumed  that  the 
nucleus  is  fixed  during  the  collision.  The  total  Hamiltonian  consists  of  the  atomic  Hamiltonian,  HA,  the  kinetic 
energy  operator  of  the  free  electron,  7^,  and  the  Coulomb  potential  V  between  the  free  electron  and  the  atom. 


H  =  HA+re  +  V, 

(2) 

N  N 

N  ^ 

]2><+Zi  1 

-IA. 

F- 

(3) 

1=1  i>j  \  '  J  \ 

i=  1  1  z| 

r  —  \ 

J  e  ~  ~  V7V+1> 

(4) 
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(5) 


Here  Z  is  the  nuclear  charge.  The  origin  of  the  coordinate  system  is  chosen  to  be  at  the  nucleus. 

The  solution  to  Eq.  (1)  are  obtained  by  first  solving  the  atomic  Schrodinger  equation. 

(HA-EAmT\L  rv)  =  0.  (6) 


In  the  range  of  incident  electron  energy  of  interest  to  hypersonic  flow  modelling,  from  threshold  to  « 100  eV, 
the  close  coupling  method  [10]  provides  the  most  reliable  solution  of  Eq.  (1). 

00 

T(r,L  TN+l)=^Jl{fm(TN+l)<&m(Ti\-  Tv)}-  (7) 

m=  0 


The  antisymmetrizer  J4  permutes  the  free  electron  with  the  bound  electrons  to  account  for  the  fact  that 
electrons  are  indistinguishable  and  must  satisfy  Fermi  statistics.  The  summation  is  over  all  possible  states  of 
the  atom,  including  the  continuum.  The  summation  in  Eq.  (7)  becomes  integration  in  the  continuum  region. 
The  excitation/ionization  cross  section  is  obtained  by  analyzing  the  asymptotic  behavior  of  the  function  fin. 
Because  there  are  an  infinite  number  of  discrete  and  continuum  target  states,  for  practical  reasons  the 
summation  is  necessarily  truncated.  The  success  of  a  close  coupling  calculation  depends  on  the  quality  of  the 
atomic  wave  functions  used  and  the  number  of  terms  included  in  the  summation.  Two  of  the  most  successful 
approaches  used  in  e-atom  collisions  are  the  convergent  close-coupling  (CCC)  method  [11,  12,  13]  and  the  R- 
matrix  with  pseudo  states  (RMPS)  [14,15].  In  the  CCC  method  expands  target  states  using  square -integrable 
functions.  The  convergence  of  such  representation  is  tested  by  successively  increasing  the  size  of  the  basis. 
While  the  CCC  method  has  demonstrated  excellent  results,  so  far  its  application  is  limited  to  atoms/ions  with 
one  or  two  electrons  outside  a  closed  shell  core,  i.e.,  the  alkali  and  alkaline  earth  atoms.  The  current  version  is 
not  applicable  to  N,  O,  and  their  ions  that  exist  in  the  Earth  entry  environment.  The  RMPS  method  is  an 
extension  of  the  R-matrix  method  [16].  The  target  functions  used  in  the  Eq.  (7)  are  constructed  using  both 
physical  orbitals  determined  from  the  solution  of  Eq.  (6)  and  pseudo  orbitals  constructed  using  a  Sturmian- 
type  basis  or  B-spline.  The  latter  approximately  represent  the  high-lying  bound  states  and  continuum  states. 
This  method  is  applicable  to  atoms/ions  of  arbitrary  structure,  but  the  accuracy  of  the  calculation  depends  on 
the  number  of  physical  target  states  used  in  the  expansion. 

Simulation  of  nonequilibrium  gas  dynamics  requires  a  complete  set  of  cross  sections  involving  all  possible 
initial  and  final  states  of  the  atom/molecule.  The  NIST  database  [18]  lists  381  levels  for  N  atom  not  counting 
the  high  Rydberg  states.  Thus  a  very  large  set  of  cross-section  data  are  needed.  As  discussed  in  the  preceding 
paragraph,  modern  quantum  mechanical  calculations  can  provide  reliable  cross  data  but  the  accurate 
calculations  are  limited  to  a  small  set  of  low-lying  states  for  a  given  atom,  or  applied  only  to  a  particular  class 
of  atoms.  Therefore  it  is  necessary  to  employ  approximate  methods.  The  most  frequently  employed 
approximations  are  classical  mechanics  or  Born  approximation  based  on  screened  hydrogenic  target  functions. 
Both  are  attractive  to  modellers  of  macroscopic  phenomena  because  analytic  formulas  are  available,  making 
the  calculation  of  large  numbers  of  cross  sections  feasible.  At  present  the  best  practice  is  to  use  experimental 
or  more  accurate  quantal  treatment  when  available,  and  supplement  the  data  set  using  either  the  classical  or 
the  hydrogenic  approximations. 
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It  is  difficult  to  predict  a  priori ,  under  what  conditions  these  approximation  will  fail.  As  an  example,  consider 
a  calculation  of  Li  atom  using  three  different  types  of  target  wave  functions,  screened  hydrogenic,  Hartree- 
Fock  (HF)  [17],  and  Multiconfiguration  Hartree-Fock  (MCHF)  [17]  functions.  In  the  HF  method,  the  electron 
sees  the  averaged  potential  from  other  electrons  and  the  atomic  wave  function  is  an  antisymmetrized  product 
of  one-electron  orbitals  (one-electron  solution  of  the  Fock  equation).  Electron  correlation  arising  from 
instantaneous  electron-electron  interaction  is  neglected.  The  MCHF  method  accounts  for  electron  correlation 
approximately  by  expanding  the  wave  function  in  terms  configuration  state  functions,  i.e,  antisymmetrized 
products  one-electron  orbitals.  The  ground  state  of  Li  has  the  electronic  configuration  ls22s  2S  and  the 
electronic  configuration  of  the  first  excited  state  is  ls22p  2P.  Since  the  Is  elections  are  more  tightly  bound  than 
the  2s  or  2p  electron,  the  screened  hydrogenic  approximation  should  apply  well  to  these  states.  Table  1 
presents  selected  properties  calculated  using  these  wave  functions  and  experimental  data  when  available.  For 
the  expectation  value  of  r  of  the  2s  orbital,  (<fhs I  >\  (jh s),  all  three  methods  are  in  good  agreement  with  each 
other.  For  the  Einstein  A  coefficient  of  the  2S  -  2P  transition,  A2s.2p,  the  HF  and  MCHF  results  are  in  good 
agreement  with  NIST  tabulation  [18],  but  A2s,2p  calculated  using  screened  hydrogenic  functions  is  ~  a  factor  of 
5  smaller.  The  electron-impact  excitation  cross  section  for  the  2S  -  2P  transition  is  calculated  at  4  eV  incident 
electron  energy  using  both  the  Born  approximation  [19]  and  the  BE  scaling  method  [20].  The  cross  section 
recommended  by  Wutte  et  al  [21]  based  on  a  compilation  of  experimental  data  is  also  presented.  The  BE 
scaling  cross  sections  calculated  using  with  MCHF  and  HF  wave  functions  are  in  good  agreement  with 
experiment,  with  the  MCHF  result  slightly  closer  to  experiment.  The  Born  cross  sections  are  approximately  a 
factor  of  3  larger.  On  the  other  hand,  the  Born  cross  section  from  the  screened  hydrogenic  calculation  is  in 
excellent  agreement  with  experiment,  whereas  the  BE  scaling  result  is  a  factor  of  3  smaller.  However,  based 
on  the  fact  that  the  screened  hydrogenic  calculation  gives  a  poor  result  for  A2S2P,  the  good  agreement  in  the 
Born  cross  section  is  due  to  fortuitous  cancelations  of  errors  in  which  the  screened  hydrogenic  approximation 
underestimates  and  the  Born  approximation  overestimates  the  cross  section. 

Table  1 :  Comparison  of  properties  of  Li  atom  using  screened  hydrogenic,  HF,  and  MCHF  wave 

functions 


Property 

Screened 

hydrogenic 

HF 

MCHF 

Experiment 

($2s  A  fas)  (au) 

4.00 

3.87 

3.86 

A2S,  2p( S  ') 

6.63E+06a 

3.75E+073 

3.87E+07 

3.72E+07b 

02S-2P  (A2)c,  Born 

36.6 

117.8 

106.6 

36. 6d 

o-2s-2P  (A2)c,  BE 

12.7 

39.7 

37.9 

36.6d 

aPhoton  wavelength  is  determined  from  MCHF  calculation. 

bNIST  data  [18], 

c4  eV  incident  electron  energy. 

‘'Critical  compilation  by  Wutte  et  al  [21], 


The  Li  atom  example  serves  to  illustrate  two  requirements  for  reliable  transition  properties,  high  quality 
atomic  wave  functions  and  a  treatment  of  collision  dynamics  that  is  applicable  to  low  energies.  In  the  present 
application,  the  MCHF  method  is  used.  Since  the  study  is  limited  to  light  atoms,  relativistic  effects  are  treated 
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by  solving  the  Breit-Pauli  Hamiltonian  [17].  The  collision  dynamics  treatment  used  in  the  study  of  ionization 
and  excitation  is  described  below. 


2.1  Electron-impact  ionization 

There  are  two  mechanisms  in  electron-impact  ionization.  The  first  is  direct  ionization  where  the  colliding 
electron  directly  detaches  a  bound  electron  from  the  atom  or  molecule. 

A+e^>A+  +  2e.  (8) 

The  second  mechanism  is  autoionization.  It  is  an  indirect  process  where  the  atom  is  first  excited  to  a  bound 
electronic  state  that  lies  in  the  continuum.  The  bound  state  is  sufficiently  energetic  that  it  automatically 
ionizes. 


A+e—>A*—>  A+  +2e.  (9) 

The  total  ionization  cross  section  is  the  sum  of  the  two  contributions. 

°1=°D1+°M-  (10) 

The  cross  term  between  direct  and  autoionization  has  been  neglected.  Different  methods  are  used  to  calculate 
direct  ionization  and  autoionization  cross  sections. 

2.1.1  Direct  ionization 

Quantum  theory  has  a  demonstrated  record  of  successful  calculations  of  the  total  ionization  cross-sections  by 
electron  impact  [22] .  The  Binary-Encounter-Dipole  (BED)  model  of  electron-impact  ionization  developed  by 
Kim  and  Rudd  [23]  has  been  used  for  a  large  variety  of  atoms  and  molecules  and  the  BED  cross-sections  are 
available  in  the  NIST  database  [24].  The  improved  Binary-Encounter-Dipole  (iBED)  model  [25]  was 
introduced  in  2001  to  improve  the  treatment  of  the  dipole  term  in  the  BED  model.  In  the  iBED  model,  the 
direct  ionization  cross  section  is  expressed  in  two  terms, 


® DI  ^ BinaryEncounter^~  ® BornDipole  (11) 

The  Binary-Encounter  cross-section  cr«,„ar)^-ncou„ter  describes  the  close  collision  between  the  free  electron  and 
bound  electron.  It  is  represented  by  a  modified  Mott  cross-section  with  the  incident  electron  energy  replaced 
by  the  average  energy  from  the  Binary-Encounter  model  [26].  The  dipole  Born  cross  section  accounts  for  the 
long  range  interaction  between  the  free  electron  and  the  target.  The  difference  between  the  iBED  and  BED 
models  lies  in  the  use  of  the  Dipole  Bom  cross-section  instead  of  the  Dipole  Bethe  cross-section.  The  former 
includes  the  effect  of  electron  shielding,  thus  improving  the  description  of  the  process  at  low  incident  electron 
energies.  The  formula  for  aDI  used  in  this  study  is  given  in  Ref.  25,  Eq.  (2.46).  All  parameters  used  in  the 
calculation  are  obtained  using  ab  initio  calculations. 
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2.1.2  Autoionization 

Autoionization  is  a  two-step  process.  The  initial  state  of  the  atom,  i,  is  first  excited  to  a  metastable  state,  m. 
The  latter  decays  either  by  ionization  or  by  radiative  decay.  Thus  the  autoionization  cross  section  is  given  by 
the  product  of  the  electron  impact  excitation  cross  section  aim  and  the  ionization  probability  of  the  metastable 
state,  Pj. 

(12) 


The  ionization  probability  is  given  by 

pi  k !,m 

m  k  +  k 

KI,m  +hR,m 


(13) 


with  kLm  and  kRm  the  ionization  and  radiative  rate  coefficients  of  state  in.  The  electron-impact  excitation  cross 
section  crim is  calculated  using  the  BE  scaling  method  [20]  discussed  in  Sec.  2.2. 


2.2  Electron-impact  excitation 

The  Born  approximation  [19]  is  the  most  commonly  used  method  for  the  calculation  of  electron-impact 
excitation  cross  sections.  It  is  based  on  a  perturbation  series  expansion  of  the  transition  matrix  element, 

Tif  =  (e,kf*rN+l®/(r1  L  rN)| V + VG <p F  +  VG (p VG ^ F + L  |e/k‘*r™<t);.(r1L  rN^.  (14) 

G^  =  (E-H0+i5)-\  (15) 

H(=HA+re.  (16) 


8  is  a  small  positive  number.  The  first  term  in  Eq.  (14)  is  the  Born  term.  The  series  converges  when  the 
interaction  potential  V  is  significantly  smaller  than  the  incident  electron  energy.  Thus  the  Born  approximation 
is  a  high  energy  approximation.  The  Bom  cross  section  is  given  by 


4  n 


V  J_ 
J  ~K4 


®/(^ iL  %) 


If 


1=  1 


®/(riL 


dKl 


(17) 


The  momentum  transfer  K  is  given  by  K  =  kt  -  kf,  with  k!  and  kf  the  momentum  vector  of  the  incident  and 
scattered  electron.  For  hydrogenic  functions,  the  Born  cross  section  is  given  in  closed  form. 

The  energy  regime  of  interest  in  modelling  hypersonic  flow  is  too  low  for  the  Bom  approximation  to  hold. 
The  Born  approximation,  while  describing  long  range  interaction  such  as  dipole  and  quadrupole  interaction 
well,  does  not  account  for  electron  exchange,  the  polarization  of  the  target  electrons  by  the  free  electron,  and 
the  distortion  of  the  free  electron  by  the  target  that  are  described  by  higher  order  terms  in  the  perturbation 
series.  The  BE  scaling  method  by  Kim  [20]  incorporates  these  effects  approximately  by  energy  scaling. 
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BE  Born 

—  (7 


E  „  +  B  +  Ey. 


(18) 


Here  Ee  is  the  free  electron  energy.  B  is  the  ionization  energy  of  the  initial  state,  i.e.,  the  binding  energy  of  the 
outermost  electron.  Ex  is  the  excitation  energy.  The  BE  scaling  method  has  been  used  successfully  in  treating 
electron  collisions  with  neutral  atoms  and  ions  [20,27],  For  ionic  target,  the  plane  wave  in  Eq.  (14)  is  replaced 
by  a  Coulomb  wave.  It  should  be  noted  that  this  method  does  not  include  resonance  effects,  that  is,  the 
enhancement  of  the  collision  cross  section  due  to  the  formation  of  transient  compound  states.  Resonances 
should  be  treated  separately.  The  BE  scaling  method  is  used  in  the  present  study  of  electron-impact  excitation 
and  autoionization. 


3.0  ELECTRON-IMPACT  IONIZATIN  OF  N  ATOM 

Ionization  from  ten  states  of  N  atom  is  studied.  Their  electronic  structures  and  ionization  potentials  are 
presented  in  Table  2.  The  first  three  states  share  the  ground  state  electronic  configuration,  ls22s22p3.  The 
upper  seven  states  has  one  electron  in  the  n-3  shell. 


Table  2:  Electronic  structure  and  ionization  potentials  of  10  states  of  N  atom  used  in  the  ionization 

calculation 


Electronic 

configuration 

Ionization  potential 

(eV) 

Electronic 

configuration 

Ionization  potential 

(eV) 

ls22s22p3  4S° 

14.534 

ls22s22p2(3P)3s  2P 

3.854 

ls22s22p3  2D° 

12.151 

ls22s22p2(3P)3p  4D° 

2.784 

ls22s22p3  2P° 

10.959 

ls22s22p2(3P)3p  4P° 

2.697 

ls22s22p2(3P)3s  4P 

4.208 

ls22s22p2(3P)3d  4F 

1.557 

ls22s22p2(‘D)3s  2D 

4.076 

ls22s22p2(3P)3d  4D 

1.518 

The  cross  section  calculations  use  atomic  wave  functions  calculated  using  the  ATSP2K  code  [28].  For  the 
metastable  2D°  and  2P°  states,  the  autoionizaton  rate  coefficients  are  taken  from  Kim  and  Desclaux  [29]. 
Results  for  the  three  lowest  states  have  been  reported  previously  [30]. 

Figure  1  compares  the  calculated  cross  section  with  the  cross-beam  experiment  of  Brook  et  al  [31].  The 
composition  of  the  N  atom  beam  used  in  the  experiment  has  been  analyzed  by  Kim  and  Desclaux  [29]  to  be 
approximately  70%  4S°  and  30%  2D°.  Figure  1  presents  the  calculated  electron-impact  ionization  cross  section 
for  such  mixture  of  N  atoms  and  compares  with  experiment.  The  reported  experimental  error  is  ±  40%  at  low 
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energies  and  +  15%  at  higher  energies.  As  seen  in  the  Figure,  the  calculated  cross  sections  agree  with 
experiment  to  within  their  error  estimate. 


Figure  1.  Ionization  cross-section  as  a  function  of  free  electron  energy  for  a 
mixture  of  N  atoms  composed  of  70%  4S°  and  30%  2D°.  Experimental  data  are 
from  Brook  et  al  [31]. 


Figure  2.  Ionization  rate  coefficients  for  10  states  of  N  atom  as  a  function  of 
free  electron  temperature. 
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The  ionization  rate  coefficients  from  10  states  of  N  are  presented  in  Figure  2.  It  is  seen  that  the  rate 
coefficients  separate  into  two  groups.  The  three  lowest  states,  being  more  tightly  bound,  have  smaller 
ionization  rates  than  the  seven  states  that  has  one  electron  in  the  n= 3  shell.  The  difference  is  particularly 
striking  at  low  electron  temperatures.  The  large  ionization  rates  of  the  upper  states  lends  to  the  possibility  that 
the  upper  states  will  first  reach  Saha  equilibrium  with  the  free  electron  before  they  reach  Boltzmann 
equilibrium  with  the  lower  states. 


T  (K)  T  (K) 

Figure  3.  Boltzmann  ionization  rate  coefficients  of  N  atom  as  a  function  of  free 
electron  temperature.  10  electronic  states  of  N  atom  are  used  in  the 
Boltzmann  calculation.  The  electronic  temperature  of  the  atom  is  (a)  7,500  K, 
(b)  10,000  K,  (c)  15,000  K  and  (d)  20,000  K.  The  contributions  from  the  3 
lowest  states  and  the  7  upper  states  are  also  presented. 


Since  most  CFD  simulations  employ  a  single  ionization  rate  coefficient  for  each  species,  Boltzmann  rate 
coefficients  calculated  using  a  Boltzmann  distribution  of  the  10  states  at  7,500,  10,000,  15,000  and  20,000  K 
electronic  temperatures  are  presented  in  Figures  3  a  -  d.  The  calculation  neglects  the  remaining  upper  states  of 
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the  atom  and  separates  the  internal  electronic  temperature  of  the  atom  from  the  free  electron  temperature.  Also 
shown  are  the  contributions  from  the  3  lowest  states  and  the  7  upper  states.  At  7,500  and  10,000  K,  the 
Boltzmann  rate  is  determined  from  the  ionization  of  the  upper  states  at  low  Te,  but  determined  by  the  lower 
state  rates  at  Te  larger  than  the  internal  electronic  temperature  of  the  atom.  At  15,000  -  20,000  K,  the 
Boltzmann  rate  is  determined  by  the  upper  state  rates.  Note  that  under  hypersonic  entry  conditions,  the  upper 
states  are  likely  not  to  be  Boltzmann  equilibrium  with  the  3  lowest  states  and  their  number  densities  are  likely 
to  be  orders  of  magnitude  smaller  than  the  Boltzmann  distribution.  Instead,  the  upper  states  are  in  Saha 
equilibrium  with  the  free  electron.  The  FIRE  II  analysis  by  Panesi  et  al.  [32]  provides  a  good  illustration  of 
this  situation.  Thus  a  single  ionization  rate  for  each  species,  as  presented  in  Fig.  3  for  N  atom,  most  likely  is 
not  a  good  approximation. 

4.0  ELECTRON  RECOMBINATION  WITH  N+  ION 

There  are  two  types  of  electron-ion  recombination  processes,  radiative  recombination  (RR)  and  dielectronic 
recombination  (DR).  In  both  cases  the  focus  of  the  present  study  is  in  the  radiance  produced  by  the 
recombination  because  both  RR  and  DR  can  be  a  significant  source  of  radiative  power. 

4.1  Radiative  recombination 

In  radiative  recombination,  the  electron  recombines  with  an  ion  and  a  photon  is  emitted. 

N'(nTs')+e—>  N(nls)  +  hv.  (19) 

It  is  the  reverse  of  photoionization.  The  photon  energy  emitted  in  RR  is  given  by 

hv=EN+-EN+Ee.  (20) 

The  process  is  illustrated  schematically  in  Figure  4.  Because  the  electron  energy  is  continuous,  the  radiation 
emitted  is  also  continuous. 


Photon 
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Electron 


energy 

'l  enercr 
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Figure  4.  Schematic  diagram  of  radiative  recombination. 
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The  radiative  power  produced  by  RR  is  given  by  [33], 

d(^p-)  =  -T=  NeUe  J  shv(7n,lv_nls(£)exp(-£)d£.  (21) 

dt  /  V^r  0 

Here  s=E/kBTe  is  the  electron  energy  scaled  by  kBTe.  In  terms  of  the  photon  wavelength  Eq.  (21)  is  rewritten 
as 


/d(hv)\  2  he 

\^rrNeVe4^T, 


sh  v<JnW_nls(£)  exp (s)X  dX. 


(22) 


B  e  o 


Thus  the  radiance  R  produced  at  wavelength  X  by  the  RR  of  A+  ions  in  the  (n  7 ’s’)  state  to  all  neutral  states  is 

R(X)  =  (2K3!2X\BTyh2c2NlonN(Jj^cjn,,,s,_n,s((:)  exp(~4  (23) 


nls 


The  total  radiance  Rmai  is  obtained  by  integrating  over  the  photon  wavelength 


Rtotal=  ]  R(A)dA.  (24) 

0 


The  RR  rate  coefficients  tabulated  in  the  literature  [34,  35]  is  related  to  the  RR  cross-section  by 


^n'l's-nls^ e)  I —  ^ e 

V 7T 


j  SC7n,,w_nls(£)™p(-£)d£. 
0 


(25) 


The  average  RR  cross-section  [33]  is  related  to  the  rate  coefficient  by 


Tn'lV-nls(Ee)~{<7n'lW-nls(Te))  =  J  SCF„,IW_n,s(s)  exp(s)ds  =  *  k»’™ 


(26) 


The  AMDPP  website  [34]  provides  RR  rate  coefficients  of  N+  recombined  to  form  neutral  N  states  up  to  n  =  8 
and  1  =  1.  Under  the  Is  coupling  scheme,  the  calculation  covers  4  ion  states  and  406  neutral  states.  Under  jj 
coupling,  there  are  6  ion  states  and  1052  neutral  states.  Temperature  covered  ranges  from  10  -  10,000,000  K. 
Nahar  and  Pradhan  [35]  also  reported  RR  data  of  N+.  They  cover  recombination  to  20  neutral  states,  with 
highest  n  =  4,  /  =  3.  Temperature  range  is  100  -  10,000  K.  Only  one  ion  state  is  covered.  In  the  present 
calculation,  we  choose  to  use  the  AMDPP  database  because  it  covers  four  ion  states  in  Is  coupling.  Two 
modifications  were  made  to  the  AMDPP  database.  First,  the  energy  levels  of  N  and  N+  in  the  AMDPP 
database  were  determined  using  approximate  quantum  mechanical  calculations.  For  most  states,  there  are 
significant  differences  between  the  AMDPP  values  versus  the  NIST  data  [18].  Hence  all  energy  levels  in  the 
AMDPP  database  are  replaced  by  the  NIST  data  when  possible.  Scaling  is  used  when  the  corresponding  data 
are  not  available.  The  second  modification  is  to  extend  the  number  of  neutral  states  that  are  formed  by  the 
recombination  of  the  ground  (3P°)  state  of  N+.  We  generate  additional  176  N  levels  from  n= 9  to  12  and  /=0  to 
7.  The  RR  rate  coefficients  for  the  additional  states  are  assumed  to  be  the  same  as  ;i=8  levels. 
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The  calculation  [30]  uses  the  parameters  Te  =  10,460  K,  Ne  -  1.19xl016  cm  3,  and  N^+  =  1.02xl016  cm  3.  It 
covers  wavelengths  from  10  nm  to  10,000  nm.  The  radiance  from  the  four  N+  states  are  weighed  by  the 
Boltzmann  factor  and  summed.  The  electronic  temperature  of  N+  is  assumed  to  be  the  same  as  the  free 
electron  temperature.  The  calculated  radiance  spectrum  can  be  divided  into  three  regions. 

(1)  120  -  10,000  nm.  In  this  region,  the  RR  spectrum  is  low-intensity,  continuous  and  rather  structureless,  a 
typical  background  spectrum. 

(2)  50  -  120  nm.  Here  the  radiance  spectrum  is  significantly  more  intense,  as  shown  in  Fig.  5(a).  It  consists  of 
a  series  of  very  sharp  peaks,  with  a  sharp,  precipitous  rise  on  the  long  wavelength  side  and  a  slower  decrease 
on  the  short  wavelength  side.  The  onset  of  the  sharp  peaks  are  at  113.09,  101.97,  88.15,  85.25  nm.  The 
maximum  radiance  in  this  region,  at  85.25  nm,  is  512.6  W  cm  3  pm  1  sr"1.  By  comparison,  the  maximum 
radiance  in  the  120-1000  nm  region  is  0.5  W  cm 3  pm  1  sr"1,  approximately  three  orders  of  magnitude  smaller. 
The  large  radiance  and  the  continuous  nature  of  the  spectrum  make  this  important  contribution  to  the  radiative 
heat  flux.  To  better  understand  the  peak  structure,  Figure  5(b)  shows  the  individual  contribution  from  the  four 
ion  states,  without  the  Boltzmann  weight  factor.  The  RR  of  the  ground  state  ion  shows  three  peaks,  from 
recombination  to  the  lowest  three  neutral  states,  2s22p3  4S°,  2D°,  and  2P°.  The  RR  of  the  lD°  ion  has  two  peaks, 
from  recombination  to  the  2s22p3  2D°  and  2P°  states  of  N.  The  RR  of  the  'S°  ion  has  one  peak  for 
recombination  to  the  2s22p3  2P°  state  of  N.  The  RR  of  the  5S°  ion  has  two  peaks  from  recombination  to  the 
2s22p3  4S°  and  2s2p4  4P  of  N.  The  sharp  onset  of  each  peak  is  characteristic  of  the  RR  process.  As  seen  in 
Figure  5(b),  the  maximum  of  RR  occurs  at  threshold  electron  energy,  after  which  the  rate  goes  down  with 
electron  energy.  This  is  reflected  in  the  decrease  of  radiance  at  the  short  wavelength  side  of  the  peak. 


Figure  5.  (a)  Radiance  spectrum  from  RR  of  N+  in  the  50-  120  nm  region,  (b) 
Radiance  spectrum  from  the  RR  of  four  lowest  N+  states. 
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(3)  10  —  50  nm.  The  radiance  spectrum  in  this  region  is  extremely  weak.  Increasing  the  number  of  ion  state  used  in  the 
calculation  may  improve  the  spectrum.  However,  the  higher-lying  ion  state  has  a  small  Boltzmann  weight  factor. 
Therefore  they  will  not  change  the  results  in  a  meaningful  way. 


4.2  Dielectronic  recombination 

In  dielectronic  recombination,  the  electron  recombines  with  an  ion  and  produces  a  metastable,  excited  state  of 
the  neutral.  The  metastable  state  can  either  emits  an  electron  and  goes  back  to  the  ion  (Auger  process)  or  emits 
a  photon  and  decays  to  a  lower  state. 

A+(n'l's')+e<rJ> A  (n  l  s  )—> A (nls)  +  hv.  (27) 

hv=  EA,-EA.  (28) 

Ee  =  EA.,-EA.  (29) 

The  emitted  photon  is  from  a  bound-bound  transition.  Thus  dielectronic  recombination  produces  a  discrete 
spectrum.  A  schematic  diagram  of  dielectronic  recombination  is  shown  in  Figure  6. 


Figure  6.  Schematic  diagram  of  dielectronic  recombination. 


The  radiance  R  at  wavelength  /,  produced  by  all  dielectronic  recombination  processes  is  given  by 

7?(T)  =  (2  K3/2A3kBTe)~1h2c2NionNeve^J£CJDRnTs'-nis(£)  exp(-ff).  (30) 

nls 

The  dielectronic  recombination  cross  section  aDRn'iw-nis  from  ion  state  ( n’l’s ’)  to  form  a  neutral  state  (nls)  is 
given  by 
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Here  the  summation  is  over  all  metastable  states  Nj  that  participates  in  the  recombination  process. 

The  dielectronic  recombination  of  N+  has  been  studied  at  Te  =  10460  K,  Ne  =  1.19xl016cm3  and  Nn+  = 
1.02xl016cm~3.  The  metastable  excited  states  N  used  in  the  calculations  are  the  series  of  states  that  converge 
to  the  2s22p2(1D)  state  of  N+:  2s22p2(*D)4s  2D,  2s22p2('D)3d  2F,  2s22p2(ID)3d  2D,  and  2s22p2(*D)3d  2P.  The 
Auger  rate  coefficients  in  the  cross  section  calculation  are  from  the  AMDPP  database  [36]  and  the  radiative 
rate  coefficients  are  from  the  NIST  database  [18]. 


Figure  7a  compares  the  total  radiance  from  DR  and  RR  in  the  98  -  112  nm  region.  The  different 
characteristics  of  the  two  recombination  processes  are  clearly  evident.  Except  for  the  sharp  threshold  due  to 
the  onset  of  the  recombination  of  the  ground  state  N+(3P°)  to  the  N(2D°)  state  at  102  nm,  the  RR  spectrum  is  a 
continuum  spectrum.  The  DR  spectrum,  on  the  other  hand,  is  composed  of  sharp  peaks  typical  of  bound- 
bound  transitions.  This  is  in  accordance  with  the  physics  of  DR  recombination  where  radiation  is  emitted  from 
the  bound,  metastable  state  of  the  neutral  atom,  N  ,  to  a  lower  bound  state.  At  the  sharp  peaks  in  the  DR 
spectrum,  the  RR  spectrum  is  lower  in  intensity  by  a  couple  orders  of  magnitude  or  more.  However,  the 
continuum  character  of  the  RR  spectrum  means  the  total  radiance  from  RR  is  significant.  Indeed,  in  the  VUV 
region,  50  -  120  nm,  the  two  contributions  are  comparable:  the  total  RR  radiance  is  6.50  W  cm  3  sr  1  whereas 
the  DR  radiance  is  7.28  W  cm"3  sr"1.  In  the  visible  and  near  infrared  region,  the  contributions  from  both 
processes  drop  significantly.  For  example,  in  the  1000  -  1500  nm  region,  the  total  radiance  from  DR  is  0.14 
W  cm 33 sr"1.  An  example  of  the  DR  spectrum  in  the  infrared  region  is  presented  in  Figure  7b.  In  this  region, 
the  RR  spectrum  is  too  weak  to  be  discernable. 
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Figure  7.  (a)  Radiance  spectrum  from  DR  and  RR  of  N+  in  the  98  -  112  nm 
region,  (b)  Radiance  spectrum  from  the  DR  of  N+  in  the  1000  -  1250  nm 
region.  RR  spectrum  is  not  discernable  in  this  region. 
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5.0  ELECTRON-IMPACT  EXCITATION  OF  N  ATOM 

Experimental  data  of  electron-impact  excitation  of  N  atom  is  sparse.  Measurements  of  excitation  cross  section 
from  the  ground  4S°  to  ls22s22p2(3P)3s  4Phave  been  reported  by  Spence  and  Burrow  [37]  and  by  Doering  and 
Goembel  [38].  Spence  and  Burrow  [37]  and  Doering  and  Goembel  [39]  reported  cross  section  measurements 
for  the  4S°  to  ls22s2p4  4P  transition,  and  Yang  and  Doering  [40]  measured  cross  sections  for  the  4S°  to  2D° 
transition.  A  number  of  R-matrix  calculations  have  been  reported  for  the  e-N  system,  but  so  far  no  RMPS 
calculations  have  been  carried  out.  Tayal  and  Beatty  [41]  reported  a  study  of  the  4S°  to  2D°  excitation  in  an  11- 
state  close  coupling  approximation  using  the  R-matrix  method.  Frost  et  al.  [42]  carried  out  a  33-term  close 
coupling  calculation  using  the  R-matrix  method  and  determined  collision  strengths  from  the  3  lowest  states  to 
21  states  with  n< 3.  They  also  reported  experimental  rate  coefficients  determined  using  spectroscopic 
measurements  in  an  arc  chamber.  Comparison  between  calculation  and  experiment  shows  agreement  for  some 
transitions  and  deviations  in  others. 

As  discussed  in  Sec.  2,  quantal  treatments  based  on  the  close  coupling  formulation  using  high-level  atomic 
wave  functions  are  not  yet  available  to  cover  high-lying  excited  states,  such  as  nitrogen  with  one  electron  in 
the  «>4  shell.  The  BE  scaling  method  is  considered  a  possible  alternative  approach  because  it  retains  the 
simplicity  of  the  Born  approximation  but  included  higher  order  effects  not  covered  in  the  Born  approximation 
by  a  scaling  method.  Figure  8  presents  cross  sections  for  the  ls22s22p3  2P°  -  ls22s2p4  2S  transition  calculated 
using  both  the  Born  and  the  BE  scaling  method.  MCHF  atomic  wave  functions  determined  using  the  ATSP2K 
code  are  used.  Since  the  2S  state  lies  above  the  first  ionization  limit,  it  has  not  been  treated  by  the  R  matrix 
method.  Experimental  data  is  also  not  available.  It  is  seen  from  Fig.  8  that  the  Bom  cross  section  is 
significantly  larger  than  the  BE  scaling  cross  section  at  low  energies  because  the  lack  of  polarization  and 
distortion  effects  in  the  Bom  treatment.  Note  also  that  the  Bom  cross  sections  in  Fig.  8  are  calculated  using 
MCHF  wave  functions.  Further  deviation  will  be  found  if  hydrogenic  wave  functions  are  used. 


Electron  energy(eV) 

Figure  8.  Electron-impact  excitation  cross  section  for  the  1s22s22p3  2P°  - 
1s22s2p4  2S  transition  of  N  atom  calculated  using  the  Born  approximation 
and  the  BE  scaling  method. 
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6.0  CONCLUSION 

In  this  study  we  have  investigated  quantum  mechanical  calculations  of  electron-collision  data  with  the 
purpose  of  generating  a  database  for  modeling  hypersonic  entry  flow.  A  major  difficulty  encountered  is  that 
simulations  of  nonequilibrium  gas  dynamics  require  a  complete  set  of  cross  sections  involving  all  possible 
initial  and  final  states  of  the  atom  whereas  modern  quantum  mechanical  calculations  can  provide  accurate 
cross  section  data  only  for  a  small  set  of  low-lying  states  of  an  atom,  or  applied  only  to  a  particular  class  of 
atoms.  Thus  a  major  goal  of  this  study  is  to  investigate  approximate  methods  that  not  only  can  extend  the 
database  beyond  what  is  currently  available  from  experiments  or  R-matrix  calculations,  but  also  can  provide 
better  quality  data  than  the  Born  approximation  or  classical  mechanics  treatment  used  in  many  current 
databases.  Electron-nitrogen  atom  collisions  are  used  as  an  example.  For  ionization,  a  combination  of  the 
iBED  method  for  direct  ionization,  and  BE  scaling  method  for  autoionization  provides  ionization  cross 
sections  of  N  atom  in  good  agreement  with  experimental  data  for  the  low-lying  states.  The  method  is  then 
used  to  calculate  the  ionization  cross  sections  for  higher-lying  states.  It  is  seen  that  the  ionization  rate 
coefficients  for  individual  states  separate  into  two  groups.  The  three  lowest  states,  being  more  tightly  bound, 
have  smaller  ionization  rates  than  the  seven  states  that  has  one  electron  in  the  n-3  shell.  The  difference  is 
particularly  striking  at  low  electron  temperature.  The  large  ionization  rates  of  the  upper  states  lend  credence  to 
the  possibility  that  the  upper  states  will  first  reach  Saha  equilibrium  with  the  free  electron  before  they  reach 
Boltzmann  equilibrium  with  the  lower  states.  It  is  likely  that,  under  such  circumstances,  a  single  ionization 
rate  for  each  species,  such  as  commonly  used  in  CFD  modeling,  may  not  be  adequate.  For  electron 
recombination,  the  radiance  from  radiative  recombination  and  dielectronic  recombination  of  N+  is 
investigated,  based  on  the  rate  coefficients  from  the  AMDPP  database.  These  two  processes  contribute 
significantly  to  the  radiative  heat  load  in  the  VUV  region,  a  region  not  easily  accessible  in  experiment.  For 
electron -impact  excitation,  the  cross  section  for  the  ls22s22p3  2P°  -  ls22s2p4  2S  transition  of  nitrogen  is  used  to 
illustrate  the  difference  between  the  Born  approximation  and  BE  scaling  method.  A  large  difference  is  found 
in  the  low  energy  region  important  in  modelling  hypersonic  flows.  This  is  ascribed  to  the  lack  of  polarization 
and  distortion  effects  in  the  Born  treatment.  Thus  the  BE  scaling  method  should  offer  a  better  description  of 
the  excitation  process. 
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8.0  NOMENCLATURE 


A 

A4 


a; 


N*N 


C 


=  neutral  atom  or  molecule 
=  ion 

=  Auger  rate  coefficient  for  the  autoionization  of  a  metastable  state 
=  Einstein  A  coefficient 

=  radiative  rate  coefficient  of  a  metastable  state  ( n  Is  )  to  the  final  state  (nls) 
=  speed  of  light 
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convergent  close-coupling  method  for  e-atom  collisions 

dielectronic  recombination 

electron 

total  energy  of  the  electron  +  atom  system 

energy  of  atom  in  (nls)  state 

energy  of  ion  in  ( n  7 ’s  ’)  state 

energy  of  atom  in  metastable  (n  Is)  state 

free  electron  energy 

ionization  energy  of  atom  in  (nls)  state 

energy  of  nitrogen  atom  in  (nls)  state 

energy  of  nitrogen  atom  in  metastable  (n  l  s')  state 

energy  of  nitrogen  ion  in  (n  7 ’s  ’)  state 

excitation  energy 

free  electron  function  associated  with  the  mh  state  of  the  atom  in  a  close  coupling  calculation 
cumulative  fraction  of  integrated  radiation  intensity 
Green’s  function 

degeneracy  factor  of  the  metastable  state 
degeneracy  factor  of  the  ion 
Planck’s  constant 
Hamiltonian 
Hamiltonian  of  the  atom 

Hamiltonian  of  the  atom  plus  electron,  without  interaction 
Hartree-Fock  method 

momentum  transfer  vector  of  the  free  electron 
Boltzmann  constant 

momentum  vector  of  the  scattered  electron 
momentum  vector  of  the  incident  electron 
electron-impact  ionization  rate  coefficient 
electron-impact  ionization  rate  coefficient  of  state  m 

radiative  recombination  rate  coefficient  of  ion  in  (n  7  ’s  ’)  state  to  neutral  in  (nls)  state 

radiative  decay  rate  coefficient  of  state  m 

angular  momentum  quantum  number  of  neutral  atom 

angular  momentum  quantum  number  of  neutral  atom  in  a  metastable  state 

angular  momentum  quantum  number  of  ion 

Multiconfiguration  Hartree-Fock  method 

principal  quantum  number  of  neutral  atom 

principal  quantum  number  of  neutral  atom  in  a  metastable  state 

principal  quantum  number  of  ion 

nitrogen  atom 

nitrogen  atom  in  a  metastable  state 

the  ith  metastable  state  of  nitrogen  atom 

nitrogen  ion 

electron  number  density 

ion  number  density 

nitrogen  ion  number  density 

ionization  probability  of  the  metastable  state  m 

quasi-steady  state  approximation 

radiance 
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total  radiance 

R-matrix  with  pseudo  states  method  for  electron-atom  collisions 

radiative  recombination 

spin  quantum  number  of  neutral  atom 

spin  quantum  number  of  neutral  atom  in  a  metastable  state 

spin  quantum  number  of  ion 

time 

electron  temperature 

transition  matrix  between  atomic  state  i  and  / 
kinetic  energy  operator  of  the  free  electron 

Coulomb  interaction  potential  between  the  free  electron  and  the  atom 

average  electron  velocity 

a  small  positive  number 

scaled  electron  energy 

atomic  wave  function 

one-electron  orbital 

photon  wavelength 

photon  frequency 

electron-impact  ionization  cross-section 
direct  electron-impact  ionization  cross-section 
autoionization  cross-section 


(^Binary Encounter—  Binary-Encounter  cross-section 
( TDipoieBom  =  Dipole  Bom  cross-section 

er„  'i ’s’ -#iis  -  radiative  recombination  cross-section  of  ion  in  (n  7  ’s  ’)  state  to  neutral  in  (nls)  state 
«7„  rs’-mis>=  average  radiative  recombination  cross-section  of  ion  in  (n  7  ’s  ’)  state  to  neutral  in  (nls)  state 
aDR  n’ly-nis  =  dielectronic  recombination  cross  section  from  ion  state  (n  7 ’s  ’)  to  neutral  state  (nls) 
i  =  Spatial  and  spin  coordinates  of  the  electron 

VF  =  Wave  function  of  the  electron  +  atom  system. 

V2  =  Laplacian  operator 
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